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ABSTRACT 



High energy emission (> tens MeV) of Gamma-Ray Bursts (GRBs) provides 
an important clue to understand the physical processes involved in GRBs, which 
may be correlated with the GRB early afterglow. A shallow decline phase has 
been well detected in about half Swift Gamma-ray Burst X-ray afterglows. The 
widely considered interpretation involves a significant energy injection and pos- 
sibly time-evolving shock parameter (s). This work we calculate the synchrotron- 
self-Compton (SSC) radiation of such an external forward shock and show that 
it could explain the well-known long term high energy (i.e., tens MeV to GeV) 
afterglow of GRB 940217. We propose that the cooperation of Swift and GLAST 
will help to reveal the nature of GRBs. 

Subject headings: gamma-rays: bursts - ISM: jets and outflows - radiation 
mechanisms: nonthermal 



1. Introduction 

Among the high energy (above tens MeV) afterglow of Gamma-ray Burst (GRB) de- 
tected so far, that of GRB 940217 is the longest and also the most energetic one. The 
sub- GeV emission lasted more than 5000 seconds and it included also a 18 GeV photon 



( iHurley et al. 19941 ) . The spectrum in the energy range 1 MeV to 18 GeV, cannot be fitted 
with a simple power law (see Fig. 3 of Hurley et al. 1994). A new spectral component in 
the energy range larger than several tens MeV is needed. This finding motivates many inter - 



esting ideas: (i) the interaction of ultra-relativistic protons with a dense cloud (IKatz 19941 ) 
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ii) synchrotron-self-C ompton (SSC) scatteri ng in early forward/reverse shocks or during 



the prompt emission (IMeszaros Rees 19941 ) a nd (hi) an e lectromagnetic cascade of TeV 



7— rays in the infrared/microwave background (IPlaga 19951 ). However, the two important 
observation facts: (a) the count rate of high energy photons is almost a constant (b) the 
typical energy of these photons is nearly unchanged, have not been satisfactorily reproduced. 

To interpret the unusual high energy afterglow of GRB 940217, the physical processes 
involved in the early GRB afterglow phase is highly needed. The successful launch of the 
Swift satellite did open a new window to reveal what happens in the early GRB afterglow 
phase. As summarized in Zhang et al. (2006) and Nousek et al. (2006), in a canonical 
Swift GRB X-ray afterglow lightcurve some interesting features are emerging, including the 
very early sharp decline (i.e., phase-I), a shallow decline of the X-ray afterglow (i.e., phase- 
II), and the energetic X-ray flares (i.e., phase- V). The interpretation and the implication of 
these features have been discussed in great detail (see Meszaros 2006; Piran & Fan 2007; 
Zhang 2007 for recent reviews). For Phase-II, which interests us here, a widely considered 
explanation is a significant energy injection dE-^/dt oc t~ q (see Zhang et al. 2006 and 
Nousek et al. 2006 and the references therein). However, the energy injection process, if 
there is, seems to be not all the story. As shown in Fan & Piran (2006), for some GRBs 
with good quality multi-wavelength afterglow data, the X-ray and optical light curves break 
chromatically and thus challenge the energy injection model (see also Panaitescu et al. 2006, 
Huang et al. . 2007). An assumption additional to the energy injection to solve such a puzzle 
is that the shock parameter (s) may be shock strength dependent (i.e., time-dependent). 

In this work, we calculate the synchrotron-self-Compton (SSC) radiation of the external 
forward shock undergoing a significant energy injection. The shock parameters, e e and e^, 
the fraction of shock energy given to the shocked electrons and magnetic field, are assumed 
to be time-dependent. We show the high energy afterglow of GRB 940217 could be given 
rise to in such a scenario. 



2. The SSC emission of the forward shock 



In this section we calculate the synchrotron-self-Compton (SSC) radiation of the external 
forward shock undergoing a significant energy injection and evolving shock parameters. 
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2.1. The SSC emission in the standard fireball model 



Synchrotron radiation. As usual, we assume that the electrons accelerated by 
the shock would follow the power law distribution dn e /d^ e oc 7~ p for 7 e > 7 m , where 
7 m is the minimum Lorentz factor of shocked electrons (Sari et al. 1998). Then follow- 
ing Yost et al. (2003) and Fan & Piran (2006), the observed typical frequency of syn- 
chrotron radiation is v m = 4.2 x lO 1 ^) 1 /^ -^B-^^h 3/2 Hz - Where z is the redshift 
of the GRB, e e and e B are the energy fraction occupied by the electrons and magnetic 
field respectively, E is the isotropic energy, t is the observer time. Here the convention 
Qx = Q/^0 X has been adopted in cgs units throughout the text, and we have taken the 
spectral index of the electron distribution p = 2.5. The observed cooling frequency is 

— 3 

v c = 4.1 x 10 16 (^ £ )~ e B _ 2 E 5 2^ 2 n~ 1 t 3 1//2 (1 + Y)~ 2 Hz, where n is the surrounding medium 
density, Y = [—1 + ^/l + Axe e /e B ]/2 is the Compton parameter, x ~ min{l, (f m /f c )^ _2 ^ 2 } 
(Sari & Esin 2001). The peak flux is F Vjfnax = 3.2(^)e^ 2 _ 2J B 52 n 1 / 2 L'^ 2 mJy. Where D 28 
is the luminosity distance in units of 10 28 cm. So in the fast cooling phase [y c < u m ) the 
light curve is (Sari et al. 1998): F v oc t^ 6 (l + F) 2 / 3 for v c > u, F v oc rV 4 (l + Y)' 1 for 
v m > v > v c , and F v oc £~( 3 p~ 2 )/ 4 (1 + Y)^ 1 for v > u m . While in the slow cooling phase 
{y c > v m ) the light curve is: F u oc t 1 / 2 for v m > v, F v oc t -3 ^ -1 )/ 4 for v c > v > u m , and 
F u oc t-( 3 P- 2 )/ 4 (l + F)- 1 for v > v c . 

Please note that in previous works the evolution of the Compton parameter Y has always 
been ignored. However, If Y 3> 1, then the effect of Y evolution should be considered. From 
the relation x = mm{l, (^m/^c)^ -2 ^ 2 } we know that, in the fast cooling phase, x — 1, then 
Y is independent of the time, while in the slow cooling phase Y should evolve with time. 

x = (is c /v m )-(P- 2 V 2 oc (1 + y)p-2 r (p-2)/2 > then we get 



y ~ 10^(^)^e:%e B %Ei^n^f 3 ^ (1) 

So we find Y oc t"(p- 2 )/ 2 ( 4 -p) , then v c oc #p-8)/ 2 ( 4 -?>) , an d so we have F v oc * VAT" 4 for 
z/ > z/ c . 

The SSC emission. The effect of SSC process on GRB afterglow emission has been 
discussed by several authors (e.g. Wei & Lu 1998, 2000; Sari & Esin 2001). The typical 
frequency of SSC emission is 

^ = H> m = 3-36 x 10 21 (1±^/ 2 

<-iWVV 4 Hz (2) 
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Where 7 m = e e {^)^T ~ 2 x 10 3 e e -i^f/i" 1 ^ 3/8 is the minimum Lorentz factor 
of shocked electrons, T is the bulk Lorentz factor, the cooling Lorentz factor 7 C ~ 2 x 
10 4 e^^ 2 E 52 3/8 n- 5 / 8 (i±^)- 1 / 2 tJ /8 (l + y)" 1 , then the cooling frequency of SSC emission is 



v f « 2 7 > c = 3.28 x 10 25 (i±^)- 3 / 2 



eZlE^n-V%^(l + Yy*Kz (3) 



The peak flux of the SSC emission is (Sari & Esin 2001) 

* \ru7TTF Vimax = 1.2 x 10~ 13 (i±^) 

i^V^^t ergcm^s^MeV- 1 (4) 

Therefore, in the fast cooling phase, the light curve is 

r f V3( 1+ y)4/3 > z/ c /c >z/, 

f„ oc J *V8(i + y)- 2 , y^>y> y{ c , (5) 

[ r (9 P -10)/8( 1 + y)-2 > I/ >^C 

and in the slow cooling phase, the light curve is 

f t, v% > v, 

F u oc I r^- 11 )/ 8 , z/ c > z/ > z/f, (6) 

Again, if considering Y evolution, then, F v oc t(9p 2 -38 P +24)/8(4- P ) for y > v w 



2.2. The SSC emission with energy injection and evolving shock parameters 

A. The SSC emission with energy injection. In the standard fireball model, the 
shock energy E is assumed to be constant. However, there are increasing evidences that 
the shock energy may increase with time during some period. One good example is the 
discovery of the " shallow decay phase" in the early X-ray light curves of many GRBs, which 
is usually attributed to the energy injection (Zhang et al. 2006; Nousek et al. 2006). Now 
we consider the case that there is significant continuous energy injection into the fireball, so 
the fireball decelerates less rapidly and the afterglow emission will show a shallow decline. 
The dynamical evolution and the synchrotron radiation signature for energy injection have 
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been discussed by many authors (Rees & Meszaros 1998; Dai & Lu 1998; Cohen & Piran 
1999; Zhang & Meszaros 2001a; Zhang et al. 2006; Nousek et al. 2006; Fan & Xu 2006). 



Here we assume the energy injection takes a form dE in Jdt oat q (jCohen fc Piran 1999 



Zhang fc Meszaros 20011 ) . then the fireball energy evolves with time as E oc t 1-9 , the fireball 
radius r oc t^ 2 ~ g ^ 4: , the minimum electron Lorentz factor j m oc T oc t~^ 2+q ^ 8 , the observed 
typical frequency of synchrotron radiation v m oc 7m-Br oc t _ ( 2+(J )/ 2 ) the cooling Lorentz 
factor 7 C oc (1 + y) _1 t( 3< ?- 2 )/ 8 ; an d the observed cooling frequency of synchrotron radiation 
u c oc ^BT oc (l + r)- 2 r (2 ~ 9)/2 , the peak flux F v>max oc iV eJ Br oc t 1 ^. The typical frequency 
of SSC emission is v 1 ^ oc 7^z/ m oc t~ 3 ( 2+,? )/ 4 ) the cooling frequency of SSC emission is 
oc 7> c oc (l + F)- 4 t( 5g - 6 )/ 4 , and the peak flux of SSC radiation is F^ c ~ \na T rF v ^ max oc 
^(6-5g)/4_ Tj s i n g these relations, we can get the light curves of synchrotron radiation and SSC 
emission with energy injection. 

B. The SSC emission with evolving shock parameters. In the standard fireball 
model, the shock parameters e e and €b are assumed to be constant. However, it is also 
possible that these quantities may vary with time. Yost et al. (2003), Fan & Piran (2006) 
and Ioka et al. (2006) have discussed the afterglow emission with e e and e# being time- 
dependent. By modeling the afterglow of several GRBs, it was found that the values of e e 
and €b of the forward shock are quite different from that of reverse shock (Fan et al. 2002; 
Zhang, Kobayashi & Meszaros 2003; Kumar & Panaitescu 2003; Wei et al. 2006). We note 
that the forward shock is ultra-relativistic, while the reverse shock is mild-relativistic, so this 
result suggests that the shock parameters may be correlated with the strength of the shock. 

Following Fan & Piran (2006), here we simply assume e e oc r _a , e# oc T~ b , and since 
T oc H 2 +")/ 8 , so e e oc t( 2 +9W8, e B oc t^^/ 8 . Then we can obtain the light curve of 
synchrotron radiation and SSC emission. 

C. The SSC emission with both the above effects. It is also possible that during 
the shock evolution, both the shock energy and the shock parameters evolve with time. Based 
on the previous analysis, we can obtain the synchrotron radiation and SSC emission light 
curves easily. Table 1 gives the temporal index a of the afterglow emission, where F u oc t~ a 
is adopted. We define a = a + c*e + a v + ay, where a corresponds to the contribution of 
the standard emission, cve represents the contribution of the energy injection, a v stands for 
the contribution of evolving shock parameters, and ay comes from the evolution of Compton 
parameter Y. For example, if we only consider energy injection, then a = «o + ole- If only 
evolving shock parameters is considered, then a = a + a v . If both the effects are considered, 
then a = a + cue + ct v . HY^§>1, then the term ay should be included. 
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3. The case of GRB940217 

GRB940217 was a very famous burst for its long-lasting high energy afterglow emission 
(Hurley et al. 1994). The high energy photons (E > 30 MeV) were recorded for about 5400 
seconds, including an 18 GeV photon ~ 4500s after the low energy emission had ended. The 
30 MeV to 30 GeV EGRET spectrum is (1.3±0.4)x 10- 8 (£/86Me\/)~ 2 - 83±0 - 64 photons cm- 2 s- 1 MeV -1 , 
excluding the 18 GeV photon. By integrating this spectrum, the ffuence at > 30 MeV is 
7 x 10~ 6 erg cm~ 2 . In addition, for this GRB there are two important observation facts: (a) 
the count rate of high energy photons is almost a constant; (b) the typical energy of these 
photons is nearly unchanged. This two facts imply that the flux should be nearly a constant 
~ 1.4 x 10~ 9 erg cm~ 2 s _1 . 

Based on the above analysis, if we assume that the energy injection occurred at time 
t ~ 500s, and we take the parameters as follows: e e ,_i ~ 0.7, 6b,-2 ~ 0.5, n ~ 1, ~ 5, 
z ~ 0.1, then at this time, the typical frequency of SSC emission is v 1 ^ ~ 20 MeV, the 
cooling frequency of SSC emission is v^p ~ 60 GeV, and the peak flux of SSC emission is 

~ 1.4 x 10~ n ergs cm _2 s _1 MeV _1 , which is quite agreement with the observation. In 
addition, we note that the observed photon energies lie between and z/^ c , and since the 
spectrum is very soft, (3 = 1.83 ± 0.64, so the index of electron distribution is p ~ 4. From 
table 1 we can find that, in the standard case, the flux would decrease with time as t" 25 / 8 , 
which is obviously inconsistent with the observation. If we consider energy injection, then 
F v oc t-( 6 + 19 9)/8 ; since < q < 1, so the flux would decay steeper than t~ 3 / 4 , which is still 
inconsistent with the observation. Then, there is only one choice - the shock parameters 
should evolve with time. If we take q ~ 0.5 (which is suggested by the recent Swift XRT 
data, see Zhang et al. 2006), a ~ 1, b ~ 0, then the flux would be nearly a constant, which 
is agreement with the observation. 

In order to investigate the SSC emission more carefully, Fan et al. (2007) have developed 
a numerical code to calculate the Compton process self-consistently for GRB high energy 
afterglow emission. Using this code, we calculate the SSC emission numerically, the result 
has been shown in Fig.l, the parameter are: the initial kinetic energy is 4 x 10 52 erg, n — 1, 
z = 0.1, eB,-2 = 0.3, p = 4, and 9j = 0.2. For 10 s < t < 5000 s, the energy injection form 
is taken to be dE/dt = 4 x 10 50 (t/10s)~°- 55 and e e = 0.06(t/5000s) a4 . At late times, the 
energy injection disappears and e e = 0.06. From Fig.l we find that the numerical results are 
consistent with the analytic estimates and can well account for the observation (both the 
light curve and the spectrum) of GRB940217. 
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4. Discussion and conclusion 

Since its discovery, the long-lasting high energy afterglow emission of GRB 940217 has 
been extensively discussed (e.g. Katz 1994; Meszaros & Rees 1994; Plaga 1995; Cheng & 
Cheng 1996; Dermer et al. 2000; Wang et al. 2001; Zhang & Meszaros 2001b; Guetta & 
Granot 2003; Pe'er & Waxman 2004). However, we note that the SSC emission of the stan- 
dard afterglow model cannot well account for the observation of GRB940217. For example, 
Dermer et al. (2000) calculated the synchrotron-self-Compton emission during the blast 
wave propagation, but their spectrum is very hard, which is inconsistent with the observa- 
tion of GRB 940217. As shown in Tab.l, the flux of SSC emission would decrease with time 
as t~ 25 / 8 for a p ~ 4 that is needed to reproduce the very steep MeV to GeV spectrum, 
even when considering the energy injection, the flux would still decay steeper than t~ 3//4 . 
Therefore, the nearly constant of the high energy flux strongly suggests that other physi- 
cal processes should be involved, such as the evolution of shock parameters and/or energy 
injection considered in this paper. Some possible energy injection processes have been pro- 
posed in Rees Sz Meszaros (1998) and Dai & Lu (1998), while the possi ble physical scenari o 



giving rise to the time-evolving shock parameters are far from clear (IPiran fc Fan 20071 ). 



The peculiar chromatic break detected in q uite a few early X-ray and optical afterglow data 



(IFan fc Piran 20061 ; iPanaitescu et al. 20061 ). however, did indicate such a possibility. 



In previous analysis, the dependence of the Compton parameter Y on time has always 
been ignored. However, from table 1 we can see that, in some cases, the influence of Y 
evolution cannot be ignored. For example, in the standard case (g = 1, a — b — 0, i.e. with- 
out energy injection and the shock parameters are constant), the light curves of synchrotron 
radiation {y > u c ) and SSC emission {y > v* c ) will be flattened by t 1 / 6 and t 1 / 3 respectively 
for p = 2.5. If p — 3, then the effect will be more prominent, the light curves of synchrotron 
radiation {y > v c ) and SSC emission {y > will be flattened by t 1 / 2 and t 1 respectively. 
So under some circumstances, the effect of Y evolution should be considered. 

GLAST will be lunched soon, and it is expected that GLAST will detect high energy 
emission (20 MeV to 300 GeV) of GRBs, which may open a new window to understand the 
physical processes occurred in GRBs. We hope that GLAST can detect more events like 
GRB940217, and this can provide important clues to explore the nature of GRBs (see Fan 
et al. 2007 for extensive discussions). 
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Table 1: The Temporal index a of afterglow emission. Here F u oc t~ a is adopted. We define 
a = ao + as + ot v + ay, where «o corresponds to the contribution of standard emission, 
represents the contribution of energy injection, a v stands for the contribution of evolving 



shock parameters, and ay comes from the evolution of Compton parameter Y. 
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Fig. 1. — The SSC radiation of the forward shock undergoing energy injection and with 
evolving shock parameters, the case of GRB 940217: the thick solid line is the light curve 
and the inserted plot is the spectrum (the times have been marked in the plot), in the energy 
range of 30 MeV - 30 GeV. For 10 s < t < 5000 s, the energy injection form is taken to be 
dE/dt = 4 x 10 50 (t/10s)-°- 55 and e e = 0.06(t/5000s) a4 . At late times, the energy injection 
disappears and e e = 0.06. Other parameters involved in the calculation are as follows: the 
initial kinetic energy is 4 x 10 52 erg, z = 0.1, eB,-2 — 0.3, p = 4, and 6j = 0.2. 



